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I. INTRODUCTION 

Liouvillc space theories such as the Mori-Zwanzig 
theory may be used to describe the evolution of a 
set of quantum mechanical observable*^ under con- 
ditions of non-equilibria. To the author's knowl- 
edge, the current formalism published in the nuclear 
magnetic resonance (NMR) literature^^ does not 
include the possibility for time-dependent Hamilto- 
nian dynamics^. The usual quantum mechanical 
master equation obtained from a perturbation ap- 
proach^ is the tool of choice for describing time- 
dependent dynamics^. These methods normally 
rely on the assumption of a short correlation time. 

The Mori equations provide a convenient frame- 
work for describing time-evolution in more general 
and not necessarily perturbative cases, either by an- 
alytical or numerical solution. Modern experiments 
increasingly consist of rapid stochastic or periodic ir- 
radiations^^ for which the timescale of the applied 
pulses may approach that of molecular motions^. It 
is therefore sensible to extend the current formal- 
ism so these modern experiments can be described 
adequately. 

We begin by summarizing the usual derivation of 
the Mori equations for time-independent Liouvillian 
dynamics. We then derive the equations in the case 
case of time-dependent Liouvillian. 



II. TIME-INDEPENDENT 
HAMILTONIANS 



dA(s) 

ds 



= -e-^TrLe^ 1 ^ 1 -, 



(2) 



where L is independent of s. Integration from to t 
leads to: 



A(f) = A(0) - f e- sL 7rLe s(1 - 7r)L d4 
Jo 



(3) 



Substitution of Eq.([T|) and A(0) = 1, and multipli- 
cation on the left by e* L yields the desired result: 

e t(l- W )L = e tL _ f* e (*-s)L^ Le s(l— )L d ^ (4) 

Jo 



B. Heisenberg Equation of Motion 

The Heisenberg equation of motion for an observ- 
able A{t): 



(It 



A(t) = LA(t) 



has the following solution 



A(t) = c tL A(0). 



(5) 



(6) 



where L is the time-independent Liouvillian, defined 
by the commutator LA = i [H, A] . 



A. Dyson Formula 

The usual derivation of the Mori equations starts 
from the Dyson formulaii&£, the derivation of which 
is repeated here for convenience. Differentiation of 
the function 



A( S ) 



= c -o'L c s(l-7r)L 



(1) 



C. Matrix Operator Notation 

Consider set of n quantum operators 
{Ii, I2, ■ ■ ■ I n } and write them as a column vector 



<4 = 



h 
h 



(7) 



with respect to s gives: 
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To this column vector corresponds a row vector: (c) The Mori product, which also depends on /? 

and H: 

y = [h h ... In] ■ (8) 

Next, an inner product is chosen which is most 
appropriate to the problem under study. The nota- 
tion {A\B) denotes the inner product of A and B. 
There are several popular inner products which may 
be used, depending on the application: 

(a) The trace projectio n 10 ! 11 , also called the 
Hilbcrt-Schmidt inner product: 



Tr 



-0H 



z 



-Ah-^ H Be u ^ H 



du. 

(11) 



Next, we define the following projection superop- 
crator, 



{A\B) tmp _ = Tr [A*B] . 



(9) 



(b) The thermal average— which depends on the 
temperature (3 = 1/kT and Hamiltonian H of 
the system: 



{A\B)t.a. = Tr 



A*B- 



-/3H 



;Z = Tr[e-e H ], (10) 



7r|=r) = (jr|s?)-(s?|sf) |sf). (12) 

In component form, this is: 



3,k 



(Sf |Sf) 



(13) 



The symbol ( 'S \ C S ) denotes the n-by-n matrix 
where exp(-/3H) jZ is the canonical statistical formed by outer product of the column vector (Sf | 
operator. and the row vector |Sf ): 



h 
h 

In 



[h h ... /»]) = 



(hlh) (h\h) 
{h\h) (hlh) 



(h\In) 
(h\h) 



(I n \h) (I n \h) ... (I n \I n ) 



(14) 



whose entries are the inner products whereas 



(Sf |Sf) 

is the inverse of that matrix. Thus, a superoperator 

commutes with (Sf|^) and (Sf |Sf) because the 
entries of those matrices are c- numbers. ) | 

When the bracket is acting on a matrix & <£ , the 
inner product is taken over all entries of the matrix 
individually, as was done in Eq.fTJJ. 

We check that this is indeed a valid projection 
operator: 



Tr 2 | J?) =((Jf|Sf) • (Sfl^)- 1 ^) • (Sf|S7)~ |Sf) 



The form of the projector is particularly simple 
if the set of operators is orthogonal. For example, 
with {/+,/_, I z } the matrix (Sf |^) is diagonal: 



1+ 
h 



(I+\I+) 
[7+ I- /,])= (J_|J_) 

(J,|Z,) 

(15) 

— > <— -l 

and the inverse matrix ( Sf | Sf ) takes the diagonal 
form: 



o 





(7_|/_) 








1 



(16) 



We note also that the entries of matrices (fTSj) and 
(|16p are c-numbers and therefore commute with any 
linear superoperator. 
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D. Derivation of Mori's Equation 

The projector method separates the Liouville 
space into two subsets: 

(a) Deterministic and slow variables, 

(b) Rapid random fluctuations. 

These two orthogonal subspaces of the Liouville 
space are decomposed using projectors as follows: 
7r is the projector onto the deterministic variables 
while 1 — 7r is the projector onto the random part. 

The usual treatment&2*i2 begins from the Dyson 
formula (Eq. [p, 



Combining the left and right hand sides together 
and making the following definitions 



JT(s)=(LJf(»)|Sf)-(Sf|Sf) , (21) 

tit =(LSf|Sf) • (Sf|Sf)~\ (22) 

we get the following set of coupled integro- 
diffcrcntial equations: 



dt 



'{t))=iVL \&(t)) + \&(t)) 

+ f J?(a)\&(t-8))d8 
Jo 



(23) 



2 tL = c tlyl 7r - )L + / e^* ^TrLe^ 1 *^ L ds. (17) which read, in component form: 
Jo 



We operate with this expression on the column 
vector (1-tt)L \&). For the left hand side of Eq. pT]) . 
e* L (l-7r)L|^), we get 



— € 

dt 



J -e iL (lJ?|<#) • 1 \%) 



= | \%{t)) - ■ (V&f 1 \&(t)) , (18) 

where ' (t)) is shorthand for e' L \ For the right 
hand side of Eq. pT)) 



,t(l-7T)L 



we have: 



eC-^Le^-^d; 



(I-tt)L^) 



e t(1 - T)L (l -ir)L\&) + 

f e (t - s)L 7rLe s(1 - 7r)L (l-7r)L|^)d.s 
Jo 

--\^(t))+ f e (t - s)L 7rL|?(s))ds 
Jo 



= |#(t))+ / e ( '- s)L (L^(s)|^) • (^|^) |Sf)ds 
Jo 

= |j?(i)) + [ (L&(s)\&)-(y\2?f 1 \y(t-s))ds, 
Jo 

(19) 

where we have made the following abbreviation: 



0{t)) =e t(1 - 7r ) L (l-7r)L 



(20) 



The last equality follows from the commutativity 
of the operator e'* -s ) L with the matrices (LjF(s)|Sf ) 

and (Sf |5f) whose entries are c- numbers. 



+ f d S V^( S )|%((- S )). (24) 

J° k 

These are the famous Mori-Zwanzig equations 
- 8 ,i2 1 i3,i4 describing the time-evolution of quantum 
operators in the presence of lattice degrees of free- 
dom. 

• The first term on the right hand side is an in- 
stantaneous oscillation term. It is a Markovian 
term, in the sense that it is evaluated at the 
same time t as the left hand side of the equa- 
tion, regardless of the past. 

• The second term, J^-(i)), is a random driving 
force because it is constructed using projectors 
(1 — 7r) onto the "noise subspace" of random 
variables, according to Eq. ([20]) . This random 
force describes the effects of the bath dynamics 
on the behavior of the spin system. 

• The last term describes relaxation and is not 
necessarily exponential; rather, it is a time- 
lagged memory term giving the effects of past 
behavior on the present. Note that the quan- 
tity (\-J?(s)\&) can be viewed as the covari- 

ance of LJ^"(s) and the initial value ^, assum- 
ing stationary (in the wide sense) random pro- 
cesses. 



E. Weak Coupling Approximation 

Using definition of adjoint operator (L^ J 2 ,^) = 
, LJ2) and assuming an inner product is chosen 
such that the Liouvillian is self-adjoint i.e. = L, 
we can write the memory function (Eq. I21[) as: 
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jT(s) = {e^-^il -7r)LSf|LSf) • (Sf|Sf) (25) 
where 



C(s) = (e'< 1 - ,r > L (l-7r)I.Sf|LSf) 

is a time-correlation function. 

Following Chorini^, we consider two asymptotic 
regimes with respect to the correlation time of the 
dynamics. This correlation function contains the su- 
pcroperator, 



7rLe a(1_,r)L (l -ir)iA 



(26) 



which remains unchanged if we insert the projector 
(1-tt): 



7rL(l- 7 r)e ;i ( 1 - 7r ) L (l-7r)L^ 
This expression is also equal to 



(27) 



ttL(1 -7r)e sL (l -tt)L 

7rL(l — 7T 



L(l-7r)LSf. (28) 
Expanding the second term in a Taylor series: 

e »(i-*) _ e sL =1 + s (l - tt)L + 1 sL . . . 

= -sttL + 0(s 2 ) (29) 

and since (1 — ir)ir = 0, expression (|28|) becomes: 

7rL(l-7r)e sL (l-7r)LSf +0(s 2 ). (30) 
Integration from to t yields 



At the other extreme where C(t) decays very 
slowly, the exponential e s ( 1_7r ) L does not vary ap- 
preciably and the approximation 



e «(l-*)L = c sL 



(33) 



holds (expand the memory function in powers of s 
and retain the zeroth order term). 

In both cases, this approximation reflects the fact 
that the random fluctuations have negligible effect 
on the time-evolution. This is called the weak- 
coupling approximation^. 

In this approximation, the memory function is 
given by: 



Jf(s) = ((l-7r)e sL (l -7r)LSf|LSH ■ (Sf |#)" 



F. Short Correlation Times 



(34) 



The integral term, 



/"daY) |Sfc(f-s)) (35) 
Jo t 



can be approximated by 



k 



(36) 



if the correlation time r c is much shorter than the 
observation time t. Such approximations are accu- 
rate, for example, when the random fluctuations in 
the lattice are very rapid (t c ~ 10~ 12 s) compared to 
the observation time t ~ 10~ 3 s. 



G. Strong Zeeman Fields 



/ ttL(1 -7r)e aL (l -7r)LSfds + 0(t 3 ). (31) 
Jo 

The term 0(t 3 ) can be neglected if the correlation 
time 

t c -1 = ((1-it)L&\L&)~ 1 

x / (e s(1 - 7r ) L (l- 7 r)L^|LSf)ds (32) 
Jo 

corresponding to C(t) is short. In this case, we may 
approximate c s ( 1 ~ 7r ' L by e sL in the expression for the 
memory function. 



In this section, we work out the strong field ap- 
proximation needed to derive the famous example of 
the NMR Bloch equations. To keep the presentation 
simple, we make use of the inner product of Eq. (|9j). 
While this approach foregoes the concept of spin 
temperature and equilibrium as is normally done in 
the standard treatments of NMR relaxation^, it does 
simplify the presentation and provides some relax- 
ation rates. We writ o 10 i 16 

F = $> (37) 

i 

for the total nuclear spin angular momentum op- 
erator. For a Hamiltonian H which consists of a 
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lattice part Ul-, and a spin part subdivided into 
Zeeman H = lu F z and small perturbation Hi 
(\\Hi\\ « pToll): 



H — Hl + Hq + Hi , 



(38) 



the oscillation term for the transverse component 
F± = F x ± iF y of Sf is 

1(F±|L|F±) {F±\[H ,F±]) , (Filli/!,^]) 



z (F±|F±) 



(*±|f±) 



(*±l*±) 



(39) 

Using the commutator [F Z ,F±] = ±F± and ne- 
glecting the second term since Hi is assumed small 
compared to Hq we are left with 



H. Bloch Equations 

To recover the classical Bloch equations, we make 
the following assumptions: 

(a) The memory function is a diagonal matrix. We 
will return to this point shortly. 

(b) The correlation time is short (see Section lll Fp 

(c) The set of operators is still {F+, F_ ,F Z }, how- 
ever, in the equations of motion, we substitute 
the deviations from thermal equilibrium 

F(t) - F(t) - (F) eq . 
We then obtain the classical Bloch equations 



(F±\F±) 

For the longitudinal component we have: 



5 \F ± (t))=±ioj \F ± (t)) 



(40) 



at 
a 



T 2 



\F z (t)) 



\F z (t)) - (F z ) eq 



at' * K ~" Ti 

where the relaxation rates are 



(44) 



^ = i^J L l^.^JI^^=0. (41) 



i (F Z \F Z ) 



(F Z \F Z ) 



Thus, the f2 matrix for the set of operators 
{F + ,F_,F Z } is 



Omitting the random driving force and using these 
approximations for fl, the Mori master equations for 
F± and F z read 



1 

% 
1 

Ti 



{LF± | e* L | (1 - n)LF ± ) 

(F±\F±) 
[LF Z I e sL I (1 - ir)LF z ) 



ds, 



(F Z \F Z ) 



ds. 



(45) 



++ 


+ _ 


VL +Z 









o" 






n-z 







-Uq 







n z _ 


U,zz 














The memory function is diagonal if there are no 
couplings. The projection supcropcrator 



n\LF) = (LF\F).(F\F) \F), (46) 
takes the following form 



a_ 

at 



d_ 

at 





r(LF+ 


LF+) 


(LF + 


LF_) 


(LF+ 






(F+ 
(LF_ 


*+) 
LF+) 


(F- 
(LF_ 


F_) 
LF_) 


(^1^) 
(LF_|LF.) 


±iuj \F ± (t)) 


(F+ 


*+) 
LF+) 


C-F z 


F-) 
LF_) 


(LF, 


LF 2 ) 


pOO 


L (F+ 


*+) 




F-) 




F,) J 









F_ 




F 2 



). (47) 



3=+,-,z 



poo 

F z {t))= / ds Y, K 'A*)\F&-*)) 



3=+,-,« 



(42) 



where ify (s) is the 3x3 matrix whose elements are 



. . (Fi I Le sL (l - tt)L I F,) , . 

M • (43) 

These equations involve integral operators with 
past information. Instead of scalar relaxation rates 
1/Tx and 1/T 2 they contain convolution integrals 
over the memory functions Kij(s). Consequently, 
spin relaxation is not necessarily exponential. 



In the case of no couplings, H = Hz + Hz, and 
off-diagonal elements such as 

(LF+|LF_) = -(F+|F_) 

vanish due to the orthogonality of the F,'s. The 
matrix then simplifies to 



HtH±) 

(F+\F+) 

o 
o 











~F+~ 


-|LF_) 
-\F-) 







F_ 





(LF Z |LF 2 ) 




F z 


{F,\F.) - 







). (48) 



When couplings are weak (H-fToH <C ||i?z||), these 
simplified equations of motion [uncoupled Bloch 
equations, Eq. (|44p ] yield a reasonable approxima- 
tion. 
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III. TIME-DEPENDENT HAMILTONIANS If we write the solution in the form 

We now show how the previous treatment can be 
extended to accomodate time-dependent Hamiltoni- A(t) = U(0,i)A(0) (51) 

ans. 



then Eq. (j4*!))) implies: 



A. Time-Evolution Operator 



The starting point is the Heisenberg equation of ~ ~dJ A ^ ~ ^M ' *)^(°) ( 52 ) 
motion for an observable A(t): 

and since this relation holds for arbitrary initial con- 

, ditions A(0), we have the general formula which de- 

— Ait) = L(t)A(t) (49) fines the time-evolution operator: 
at 

where L(t) is a time-dependent Liouvillian defined 

* d -^A = L(t)U M . (53) 

L(t)A(t) = i[H(t), A(t)]. (50) The solution of this equation is given by: 
I 



u (°> t) = f] f d n r dr 2 ... r * dr n L( n )L(T 2 ) . . . L(r n ) (54) 

n=0 J ° J ° J ° 



I 

where the Liouvillians are left time ordered. Conse- is the solution of Eq. (|53[) can be seen by differenti- 
qucntly, Eq. (|49[) is a causal relation. That Eq. (fS"4")) ating Eq. (f5"4")) with respect to t: 



^ = L(i)VNri dr 2 ... dr n L(ri)L(r 2 )...L(r n ). (55) 

dt ^J Jo Jo 



I 

A shorthand notation for U(0, t) uses the left Dyson For general time evolution, we have: 
time-ordering operator T : 



U(0,t)=Texp(J^L(T)dT S \ (56) 



J 



U (*> *) = £ / d Tl [ 1 dr 2 ... [ " 1 dr„L(r 1 )L(r 2 ) . . . L(r„) (57) 



which corresponds to 



U(s,t) = T exp (J L(r)dr^ . (58) 
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B. Dyson Formula for Time-Dependent 
Liouvillian 



Now we prove the following identity, which holds, 
for time-dependent Liouvillian. 



Texp( / drL(r) J = Tcxp ( / dr(l - tt)L(t) J + f Texp 



/ \J0 J Jo 



(J drL(r)^ 7rL(s)Texp (J dr(l - tt)L(t)^ ds, 



where the right time ordering operator is defined by to the right 
the time arguments increasing as we go from the left 



Texpi^J drL(r)J = J^J d n J dr 2 . . . J dr„L(r„)L(r„_ 1 ) . . . L(n). (60) 



We first rewrite Eq. |59|) as 



-TexpQ^ rfr(l -tt)L(t)^ =-Texp^ rfrL(r) 

TexpQ| drL(r)^ ttL(s)T exp ^ dr(l - tt)L(t)^ ds, (61) 

I 

and show that each side of the equality obeys the While for the right hand side: 
same differential equation for the same initial condi- 
tions. For the left hand side, we have: 

jLHS = LHS-{l-n)\-(t). (62) 
I 



RHS = - T cxp ^ drL(r)^ L(t) + T cxp (J drL(r)^ 7rL(t) 

+ ^ TexpQf drL(r)^ 7rL(s)Texp ^ dr(l - tt)L(t)^ ds(l - 7r)L(t) 
=RHS- (1 -7r)L(f) (63) 

I 



These identical differential equations, combined proof of the Dyson formula for time-dependent Li- 
with the initial condition at t = 0, complete the 
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ouvillian. This identity is due to Holian and Evans 

17,18 



C. Derivation of the Mori Formulae Using 
Time-Dependent Dyson Formula 

Using Eq. we may derive the Mori equations 
in a manner analogous to the time- independent case. 
The analogue of Eq. (fl8|) is the following: let the left 
hand side of Eq. ||Sj5J) act on (1 - 7r)L(t) |^), 



dt 1 
where 



&i(t))-(L(t)&\&)-{&\&) |Sf f (t)) (64) 



Sft(t) = T 



exp 



drL(r) 



(65) 



is the anti-causal propagation of <& . Next, we let the 
right hand side of Eq. J5SJ) act on (1 - 7r)L(t) |#), 
and obtain an expression analogous to Eq. fT9]k 



These are the Mori equations for time-dependent 
Liouvillian. It is noteworthy that the memory kernel 
J(fj k {s,t) no longer depends on the time difference 
t — s, but also depends on the origin of time (s). 
This is a consequence of the non-commutativity of 
the Hamiltonian at different points in time. 

Another point of interest is that the dynamics are 
given for the anti-causal operators. An adjoint equa- 
tion of motion can also be derived which describes 
causal evolution, but requires adjoint definitions of 
the propagator and the expression upon which the 
Dyson equation is applied to. The weak coupling 
approximation of the previous sections remains ap- 
plicable, if it is needed. 

We also note that the Liouvillian which appears 
under the integral sign of Eq. (|7ip cannot be replaced 
by an average Liouvillian (see Eqs.l67land[68l). Such 
a substitution is not generally justified. However, 
it is possible to replace the exponential function of 
Eq. ([57|) by a Magnus expansion. 



D. Change of Reference Frames 



1^(0,*))+/ (L(s)^ t (s,i)|^)-(^|^) |^ t (s))ds 

JO 

(66) 

where 



&(8,t)) = Texp[ I dr(l - tt)L(t) I (l-Tr)L(i) \<3) 



(67) 

Combining these expressions and making the follow- 
ing definitions 



Jtri(s,t)=(L(s)&(s,t)\9)-(9\9) , (68) 

ift(t) =(L(0\&) ■ (^l&f 1 , (69) 

we get the following set of coupled integro- 
differential equations 



dt 



r(t))=in(t)|sft( t )) + |^t( 0>t )) 

Jo 



(70) 



which read in component form 



+ j\sY,^ k MWl{s)). (71) 



For convenience, we may change frames of refer- 
ence at will. This is done using time-ordered trans- 
formations, as is normally done when solving the Li- 
ouville von Neumann equation. We summarize these 
well-known transformations in the notation of this 
article. 



1. Interaction Representation 

If the time-dependent Liouvillian decomposes into 
the form L(t) = \- (t) + Ly(t), the Heisenberg equa- 
tion of motion is 



dA(t) 
dt 



[L (t) + L v (*)]A(t). 



(72) 



We define the transformed operators: 



A=Texp - / drL (r) A, (73) 



L=Texp - / drLo(r) L, (74) 



where T indicates a time-ordering from left to right, 
with the most recent Liouvillian to the right. These 
transformations arc anti- causal. Then, 
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dA d-> 
— — =— 1 exp 
dt dt 



= — T exp 
=L v (t)A. 



■J dTL (T)^jA 

■J rfrLo(r)^ L (t)A+ T exp (- J drLo(r)) [L (t)A + L v (t)A(t)] 



(75) 



Thus, in a rotating frame of reference, a Liouville 
equation 



The solution to this equation is, 



dA 
~dt 



L V A 



(76) 



holds, where Ly is the perturbation part of the Li- 
ouvillian expressed in the rotating frame. 



( J L(t')dt' 



where 



(80) 



2. Toggling Frame Transformation 

In Equation (|75p . we may further decompose the 
interaction representation Liouvillian 



L 2 = Texp - / Li(t')^' L 



(81) 



Li + L 2 



(77) 



and 



where Li is the part corresponding to the RF field 
and L2 is the remaining "internal" part of the Liou- 
villian (e.g. spin-spin interactions). 

Consider a transformation to the so-called "tog- 
gling frame" , 



= T exp ^— J drL (T 
IV. CONCLUSION 



(82) 



A = T exp 



Li (t')dt' A 



(78) 



where T is the Dyson time-ordering operator—, 
which places the most recent events to the right. 
Then, 



dA(t) 

dt 



= L 2 i(i). 



(79) 



We have extended the Mori theory to the case 
of time-dependent Hamiltonians. To the author's 
knowledge, the previously published treatments of 
Mori relaxation theory in the NMR literature did 
not include the case of time-dependent Liouvillian. 
This makes it possible to treat general cases of cw or 
pulsed irradiations in NMR and/or coherent optics, 
for arbitrary timescales of the molecular motions rel- 
ative to the period of irradiation. Also, we have out- 
lined a way of introducing coherent averaging effects 
into the Mori equations. 



* URL: http://waugh.cchem.berkeley.edu 
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